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Abstract 

We consider the hungry Volterra hierarchy from the view point of the multi boson KP 
hierarchy. We construct the hungry Volterra equation as the Backlund transformations 
(BT) which are not the ordinary ones. We call them "fractional " BT. We also study the 
relations between the (discrete time) hungry Volterra equation and two matrix models. 
From this point of view we study the reduction from (discrete time ) 2d Toda lattice to 
the (discrete time ) hungry Volterra equation. 
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1 Introduction 



The W algebra was introduced by Zamolodchikov as an additional symmetry of the con- 
formal field theory. [|I|] The classical Wn algebra can be constructed from the free fields rj 
for i = 1,2, ■ ■ ■ , N satisfying the Poisson algebra 

{r,{x),rj{y)} = {-6,,, + ^)6{x-y). (1.1) 

The Poisson structure is related A^-reduced KP hierarchy. The Poisson map from the free 
fields to the pseudo differential operator L is defined as 

L = {d + n{x)){d + r2{x)) ■■■{d + TNix)) (1.2) 
with d = d/dx, and the constraint 

N 

Y.r^{x) = 0. (1.3) 

i=l 

Recently deformations of the W algebra received much attention. One is the q deforma- 
tion of the W algebra. Another deformation of the W algebra is lattice W algebra, which 
may be related to the sIn Toda theory. The most famous example is the Lotka-Volterra 
model. It is well known that the Lotka-Volterra model reduces to the Korteweg-de Vries 
(KdV) equation in the continuous limit, and that it can be formulated in terms of lattice 
Virasoro algebra. 

In this paper we shall study the hungry Volterra equation (or " Bogoyavlensky lattice" ) 
which is known as an extended Volterra model. |Q It is pointed out that the hungry Volterra 
equation is a fundamental integrable system of the lattice algebra. We define the 
fc-hungry Volterra equation as 

dV ''■"^ ''"^ 

^ = V4T.yn+^-J:Vn-.). (1-4) 
•i^l i=l i=l 

In the case k = 1 ( |1.4|) becomes the Lotka-Volterra equation. 

We consider the hungry Volterra equation as Backlund transformations (BT) of the 
multi boson KP equation which is related to the sl{k + l, k) algebra. We call the BT 
"fractional BT" . If we repeat the fractional BT k times we can obtain the usual BT. To 
study relations between the 2d Toda lattice and the hungry Volterra equation we consider 
the hungry Volterra equation in two matrix models. The most simple model becomes 
a matrix model for the bi-colored random triangulation. Furthermore we consider the 
discrete time hungry Volterra equation in two matrix models with Penner type potential. 
It becomes a simple relation of the partition functions in bilinear form. 

This paper is organized as follows. In section 2 we consider the eigenvalue problems 
of the 3-hungry Volterra equation. In section 3 we define the fractional BT and it is 
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equivalent to the fc-hungry Volterra equation. In section 4 we consider the reduction from 
the 3-hungry Volterra equation to the classical W3 and Virasoro algebra. In section 5 we 
obtain the hungry Volterra equation and the discrete time hungry Volterra equation in 
two matrix models. The last section is devoted to the concluding remarks. 



2 Hungry Volterra equation 

We consider the spectral problems, 



A^/^vE,^ = (2-1) 



3 "3 



From (Q.l\) we can obtain 



A^^„ = ^„,+i + ao(n)^„ + aiH^„„i + a2H^„_2, 
A^„+| = ^,+|+aoH^„+2 +ai(n)^„_i +a2(n)^„_4, (2.2) 

where 

ao{n) = Un + Vn + Wn, tti (n) = f „U„_i + M„_iW„ + U7„t;„,_i , a2{n) = WnVn-lUn-2, 
do{n) = Wn+l +Un + Vn, (n) = UnWn + WnVn + t^„M„-l, a2(n) = t;„M„_iW„„i, 

do(n) = Vn+l + Wn+l + Un-, ai{n) = Wn+lVn + VnUn + Mn^n, a2{n) = M^Wnf „_i(2.3) 



Using (pT2|) we can obtain the first lattice fiow d = 



a^.+l = ^„+i + ao(n)§„,2. (2.4) 
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From ( p.2|) and ( p.4|) we can obtain the consistency conditions, 

dao(n) 



da; 



ai(n + 1) — ai(n) 



dai{n) dao{n - 1) 

- ai(n-l) + 



dx dx 

da2(n) / / N / NN / X 

= a2(ao(r2)-ao(n-2)). (2.5) 

In the same way we can obtain same equations about a{n) and a{n). These are the 
generalized To da equations. 
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On the other hands, using the consistency condition of (|2.1| ) and ( |2.4| ) the equations 
of motion read 

— — = - w„ + + 

dx 

df„ , . 



dx 

dWn 



- f„_i + M„ + (2.6) 



dx 

These are 3-hungry Volterra equations. ( |2.1| ) and (|2.4|) can be cast into the form 



where 



A^, = L,^„, A§„+i =L,+ i§„+i, A^„+| = L„+2^„+2, (2.7) 



-f/n = {d - Vn-\- Un-\){d - Vn-\- Un-\- W^~^{d - Un-\- 

{d - - - Vn)'^{d -Wn- Vn), 

Ln+\ = {d - Un^l- Wn){d - Un-l- Wn- VnY^{d - Wn- Vn) 

{d - Vn - Wn - UnY^{d - Un - Vn), 
I/n+l = {d-Wn-Vn){d-Vn-Wn-Un)~^{d-Un-Vn) 

{d-Un-Vn- Wn+l)'^{d - M„ - Wn+l). (2.8) 



This is the form of the 4-boson KP hierarchy. Using ( p.2| ) and ( p.4| ) we can obtain the 
Backlund Transformations 

Ln+i = {d - ao{n))Ln+i{d - ao)'-^, 
^n+i+i = (d - do{n))L^^{d - do{n))-\ 

Ln+i+i = (d - ao{n))L^^^{d - aoin))-\ (2.9) 

3 Fractional Backlund Transformations 

We consider multi boson KP hierarchy: 

L(^) = id- q,)id - q,)-\d - q2)id - h)'^ ■ ■ ■ {d - qk-i)-\d - q^), (3.1) 
with a Dirac constraint 

k k—1 

E^.-E^~^ = o- (3-2) 

j=l 1=1 

We construct "fractional" Backlund transformations (BT), 

(fe) ^ (fc) i^(fc) (fc) 



where 



rpik) 



Here we realize that the normal lattice jump n — n + 1 can be given a meaning 
BT, 



(3.4) 
of the 



r(fc) 
n+l+i 



and 



n+l+i 



(,) l^(fc) (fc) 



Using (0|) and (|3.6|) we can obtain a relation 



T(^) rp(K) rp(K) rp 



(fe) 



To satisfy ( p.2|) using (|3^ ) we can set 



9i 



Using ( ^.lUI ) we can rewrite ( |3.8| ) 



(f -P„,i)-'(s 



Then ( p.ll| ) becomes an equation 



dx ^1 ^- 

Using the fractional BT recursively we can obtain 



(3.5) 
(3.6) 

(3.7) 



From ( p.3|) we get 

r:?, = (a-«r*')(a-?r*V' 

To obtain the last equation we use relations 

{d-q'r'^') = {d-qt^:^^) = ... = {d-qt''-h 



From (|3.2|) we get 



(n+l+{) ^ 
i=l 



(3.8) 

(3.9) 
(3.10) 
(3.11) 

(3.12) 

k 

(3.13) 

(3.14) 
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If we substitute ( p.l4| ) into ( |3.12| ), we can obtain 



k-l 



dp. 



(3.15) 



(|3.15|) ia nothing but the /c-hungry Volterra equation. Using (|3.7] ) we can obtain the 
normal lattice BT 



n+ 



In the case A; = 3 it corresponds to (p.9|) . 



4 algebra and Torus 

We can rewrite the Lax operator (2^) using (2^) or directly from (|2.4|) , 

Ln = {d-Un-Vn-Wn)~^{d-Un+l-Vn+l-Wn+l)~'^ 
(d - Un+1 - Vn+l)id -Un- Wn+l){d - - Wn), 
{d-Un-Vn- Wn+lY^{d - U^+l - V^+l - W„+2)"^ 
{d - Un+l - Wn+2){d - Vn+l - W„+l)(5 - M„ - Vr.' 



'"n '-'nil 



-1 



{d - Vn+2 - Wn+2){d - M„+l - Vn+l){d - M„, - 

We constraint this system under the two-periodic condition 

where /j, is an arbitrary constant. Furthermore we set the total energy, 

Hl = - '^{Un + Vn+ Wn) = 0. 
n 

Using (^?T|) we can obtain 



/^^n L/n 



^(2) 



2 . 

"3 



where 



/-(2) 
^« 

r(2) 



{d - Un - Vn - Wn){d - Un+1 - Vn+l - W„+l), 
{d - Un+1 - Vn+l){d -Un- Wn+l){d - f „ - Wn), 
{d - Un -Vn- Wn+l){d - Un+1 " Vn+l " Wn), 



(3.16) 



(4.1) 



(4.2) 



(4.3) 



(4.4) 



r(l) 


— (f) 


Un+l 


— Wn)[0 — Vn+1 — Wn+1 )[0 ~ Un — Vn 




= {d- 




Vn+1 - Wn+l){d - Un+l - Vn - 


f(l) 


= (5- 


Vn - 


Wn){d - Un+l - Vn+l){d ~ Un - Wn+1 



(4.5) 

Then the Lax operators spht into two parts. The two parts are the 3 and 2 reduced type 
respectively for the constraint ( [4.31 ). It is a torus. But the equations of motion are still 
3-hungry Volterra equations. 
Furthermore we set 

^„ = ^„+i = ^, = ^„^4 = ^, =^^^1 = ^, (4.6) 

then the eigenvalue problems become 

= lW^, = lW^, A^^ = L^^, (4.7) 



where 



= (^Q — u — v){d — u — w){d — V — w), 
L^^) = (^Q — u — w){d — V — w){d — u — v), 

= {d-v-w){d-u-v){d-u-w). {U 



From ([4. 81) we can obtain the classical algebra. 



On the other hand if we set 

^ = ^ = 0, (4.9) 

the spectral problem becomes 

A^^„ = 4^)^,, (4.10) 

where 

Li^) = {^d-u„){d-Un+i) (4.11) 
( ^.11[ ) is nothing but the classical Virasoro algebra. 



5 Two Matrix Model and Hungry Volterra equation 
5.1 Hungry Volterra equation 

We consider the computation of the following integral over two Hermitian matrices U and 
V of the size N x N (see for a review) 

Z{tsk,hk,N) = J dt/dKe-^^(*'*), (5.12) 
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with the potential 

mi m,2 

V{U, V, tsk, tsk) = E i^kU^' + E ^3kV^' + 9UV, (5.13) 

k=l k=l 

and with respect to the Haar measure over Hermitian matrices H 

N 

dH = CNl[dHii n dReHijdlmHij. (5.14) 

1=1 l<i<j<N 

The parameters ^3^, tsfc, g are real numbers. We reduce the integral to an integral over the 
eigenvalues of U and V, denoted by Ui and Vi respectively. In the well known method we 
can obtain the reduced integral 

Z{t,i3k;N) = J dudvA{u)A{v)e-^^^'''^'''^''''^\ (5.15) 

where A(u) = I\i<^j{ui — Uj) and A{v) = l\i^j{vi — Vj) denote the Vandermonde determi- 
nant of matrix u and v where u = diag{ui, ■ ■ ■ ,un) and v = diag{vi, ■ ■ ■ , vn). Here even 
if we rotate the eigen values 

Ui e'^Ui, Vi C^Vi, (5.16) 

(|5.15|) does not change from the form of the potential (|5.13|) . As in the standard orthogonal 
polynomial technique, we introduce two sets of polynomials Pn = x^+lower degree, and 
Pn{y) = y^+lower degree , for n = 0, 1, 2, ■ ■ -, which are orthogonal with respect to the 
one dimensional measure inherited from ( p.l5| ), namely 

{Pn,Pm) = I dxdl/e-^(*3'='*^'=)p„(x)p„(y) = hJra,n. (5.17) 

Using the multi-linearity of the determinants, we may rewrite 

A(m) = det[M^"^]i<ij<7v = det[pj^i{ui)]i<ij<N, 

A(f) = det[vl~\<ij<N = det[pj^i{vi)]i<ij<N. (5.18) 

Using the orthogonality relations, the partition function is finally can be written as 

N~l 

Zihk, hk] N) = const, n hi. (5.19) 

Using the fact ( |5.16| ) we can obtain a relation 

Pn{xe^)pn{ye'^) =pn{x)pn{y). (5.20) 
The multiplication by x and y can be represented by the Jacobi matrixes Q and Q: 

n+l 71+1 
XPn{x) = E Qn,lPm{x), yp„{y) = E Pm{y)Qm,n, (5.21) 
m=0 m=0 
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where Qni and Qni are the matrix elements of Q and Q. 

From the definition of orthogonal polynomials it follows that 

Qn,n+1 = 1, Qn,m = 0, TTl > U + 2, Qn+l.n = 1, Qm,n = 0, TTl > U + 2, (5.22) 

Define a wave function 

^nihk, isk, x) = p„(x)e^(*3fe,t3.,-). (5.23) 
Here we introduce the matrix Q 

Qn,m i^HQH )ra,m; -^^n,m ^n^nm- (5.24) 

In what follows it will be convenient to define an explicit parameterization of matrices 
Q and Q. We choose the following parameterization 

Qn,n+1 = 1, Qn,n~k = fk{n), = 0, 1, ■ ■ ■ , 3m2 - 1, 

Qn,m = 0, for m — n > 2, and n — m > 3m2 — 1, (5.25) 

Qn,n-1 = Rn, Qn,n+k = gk{n)R~\^ ■ ■ ■ -R„+fc, /c = 0, 1, ■ ■ ■ , 3mi — 1, 

Qn,m = 0, for n — m>2, and m — n > 3m2 — 1, (5.26) 
where Rn = hn/hn-i- To satisfy the condition ( |5.2U| ) we can obtain 

/3,(^) = 0, /3,+i(n) = 0, g3,in) = 0, g^^^,{n)=0, j = 0,l,2,---. (5.27) 
We can rewrite ( |5.21| ) 



m2 



XPn = Pn+1 + f3j-l{n)Pn+l-3j, (5.28a) 

h h 

= 7 Pn-l + 2^93j-l[n)T Pn~l+3j- (5.28b) 



From ( 5.23 ) we can obtain 



x$„ = g$„, — <i>„ = q;<i>„, ^<fn = -gi$n, (5.29) 

where the subscripts "+" denotes upper triangular plus diagonal parts of the matrix and 
"-" denotes lower triangular part. 
For the first one we can obtain 

d 

— $„ = ^n+s + So(n)$„, 

OT3 

d h„ , , 

= -r^$.-3, 5.30 
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where 

soin) = h{n) + h{n + l) + h{n + 2). (5.31) 

Hereafter we consider the special case 1112 = 1. In this case the recursion relation ( ^.28a]) 
becomes 

X'^^n = + So{n)^n + Si{n)^n~3 + S2(n)<I)„_6, (5.32) 

where 

si{n) = /2(n+l)/2(n-l)+/2(n-l)/2(n)+/2(n)/2(n-2), s^in) = f2{n)f2{n-2)Mn-4). 

(5.33) 

(|5.30|) and ( |5.32| ) are nothing but the eigenvalue problem ( p.4|) and ( |2.1| ) . Then consistency 



of ( |5.3CI| ) and ( |5.32| ) becomes the hungry Vol- Terra equation 

= /2(n)(/2(n + 2) + Mn + 1) - Un - 1) - ^(n - 2)). (5.34) 

To the flows t^k we can obtain the higher order hierarchy of ( p.34|) . 
On the other hand from ( [5.30[ ) and ( |5.32| ) 



= 7^„,(so(^^) - So(n - 3)), ^^^M = 7^„+3 - 7^„, (5.35) 

where TZn = hn/hns- (|5.35|) is nothing but the 2-dimensional Toda equation. 

We consider the most simple case mi = 1712 = 1. We set g = N and ^3 = ^3 = —tN. If 
we think of t as a small parameter, we may at least formally expand the integral Z in power 
of t. This expansion is expressible as a sum over Feynmann graphs, made of double lines 
oriented in opposite directions and carrying a matrix index i G {1, 2, ■ ■ ■ , A^}, involving two 
types of three-valent vertices, connected by one type of propagator. (< UV > propagator, 
weighted by 1/A^) As usual in matrix integrals, each graph receives a contribution per 
oriented loops and tA^ per vertices. We finally obtain the expansion 

fV(r) ATxir) 

F{t, N) = log Zit, N) = Y: ^TTrFTr, (5-36) 

r |Aut(r)| 

where V^(r) denotes the number of edges, x(r) its Euler characteristic, |Aut(r)| denotes 
the order of the automorphism group of bi-colored connected graphs F. The bi-colored 
graphs mean that two adjacent vertices have different colors. We can consider ( |5.36| ) as 
an expansion over dual graphs, with now black and white alternating faces. Then we can 
finally view log Z as the partition function for black and white colored triangurations. 

5.2 Discrete time hungry Volterra equation 

In this subsection we consider following two matrix model 

Z{hk, h, l) = j df/dl^e-^^^(*'*'') (5.37) 
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with the potential 



mi 



V{U, V, hk, is, /) = E t^kU'" + hV^ + gUV - 31 log U 



(5.38) 



k=l 



We added the "log type potential" . |jTT||T^ Note that / = 0, 1,2, ■ ■ -. In the special case 
/ = 0, ( p.37| ) becomes the matrix model studied in the previous subsection. The orthogonal 
polynomials satisfy the orthogonal relations 



{Pn,hPm,l)= / dxdye ^^^^'"^^'"^^Pn,l{^)PmAy) = Kl^m,n- 



(5.39) 



The partition function (tau function) can be written 



N-l 



Z{t3,tsk,l;N) = const. Y[ Ki- 



(5.40) 



1=0 



We can obtain the recursion relations of orthogonal polynomials in the same way 

XPn,l = Pn+l,l + f2{n, l)Pn-2,l, 



yPn,l = T —Pn-l,l + 2_^gzj--l[n,l)- p„_i+3jy. 



(5.41a) 
(5.41b) 



Notice that orthogonal polynomials pn,i and pn,i and /(n, I), g{n, I) and hn,i are the func- 
tions of n and /. We rewrite /2 using the partition function ( p.40| ) 



f2{n,l) = Vn,l 



^n+3,l^n-2,l 



(5.42) 



(5.43) 



Zn+l,lZn,l 

Using ( ^.39|) and the condition e~^^*3fe:*3fe:^+i) = r^3^-vit;ik,hk,i) -^q gg^ ^ relation 

X Pn,l+1 — Pn+3,l + ^n,lPn,h ^n,l — —7 — 1^ 7^ • 

iln,l ^n-2I+l^n+l,l 

As a compatibility condition of ( |5.41a| ) and ( p.43| ) we can obtain 

In+l,l + VnA+l = Vn+3,l + ^nd, In,lVn,l = K.,«+l-^n-2,« ■ (5.44) 

Here we change dependent variables V(n,i) ^{n,i)- Notice that Vn^i = /2(n, /) are the 
dependent variables of the hungry Volterra equation ( p.34|) . We define 



Vn,l = Vn,lil + V;,«-l)(l + Vn,l-2)- 



(5.45) 



Using Vn,i we can rewrite 



[l + Vn,l){l + Vn,l-l)il + Vn,l-2)- 



(5.46) 



11 



The equation of motion ( |5.44| ) becomes 



Kl+l ^ {l + Vn-l,l){l + Vn-2,l) 

Vn,i (i + v;-i,m)(i + v;-2,m)' 



( |5.47| ) is the discrete hungry Voltrra equation. [0-[|T^ Using the partition function ( p.40[ ) 



we can rewrite 



t7 Zn+3,lZn-2,l+l , , f> Z-a iZ-a+l l+l /r Ai 

Vn,l = — , 1 + Vn^l = . (5.4J 



(|5.47|) is cast into the following bilinear equation 

Zn-2,l+lZn+3,l + -^n+l.i^n.i+l = ^n,«^n+l,/+l • (5.49) 

6 Concluding Remarks 

In this paper we consider the fc-hungry Volterra equation as Backlund transformations 
(BT) of the multi boson KP equation which is related to the sl{k + 1, k) algebra. We call 
the BT "fractional BT". If we repeat the fractional BT k times we can obtain the usual 
BT which is the Toda lattice. 

To study the relations between the 2d Toda lattice and the hungry Volterra equation 
we consider the hungry Volterra equation in the two matrix model. If we select the time 
tsk and ts from the flows of 2d Toda lattice t = (ti,t2, ■ ■ ■) and t = (^1,^2, ■ ■ ■), we can 
obtain the 3- hungry Volterra hierarchy which is related to the lattice W^- The most simple 
case is the matrix model for the bi-colored random triangulation. It is easy to construct 
the general fc-hungry Volterra equation in 2d Toda lattice. Using the duality of t and t 
we can obtain two fc-hungry Volterra hierarchies in 2d Toda lattice. 

Furthermore we consider the discrete time hungry Volterra equation in the two matrix 
model with the Penner type potential. It is well known that the partition function of 
the one (multi) matrix model with the Penner type potential satisfies the discrete time 
Toda equation. [jlO| fl^ [0 We select the flows / and ^3 where / is the multiple of 3 of 



discrete time lattice which belongs to the 2d discrete time Toda equation. If we select 
the multiple of k of discrete time lattice and tk, we can obtain the discrete time fc-hungry 
Volterra equation. Then we can find two discrete time /c-hungry Volterra equations in 2d 
discrete time Toda equation. 
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